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Abstract 

Many examples of exactly solvable birth and death processes, a typical stationary 
Markov chain, are presented together with the explicit expressions of the transition 
probabilities. They are derived by similarity transforming exactly solvable 'matrix' 
quantum mechanics, which is recently proposed by Odake and the author. The (q- 
)Askey-scheme of hypergeometric orthogonal polynomials of a discrete variable and 
their dual polynomials play a central role. The most generic solvable birth/death rates 
are rational functions of q x (x being the population) corresponding to the g-Racah 
polynomial. 

1 Introduction 

The Brownian motion, a typical stationary Markov process with a continuous state space, 
is known to be described well by the Fokker-Planck equation [TJ [2]. A birth and death 
process, on the other hand, being a typical stationary Markov chain with a set of non- 
negative integ state space [3], can be naturally considered as a discretisation of a 
one-dimensional Fokker-Planck equation. Although birth and death processes have a wide 
range of applications [21 [3], demography, queueing theory, inventory models and chemical 
dynamics, we will focus on their mathematical aspect, i.e., the exact solvability. In this paper 
we present 18 exactly solvable birth and death processes based on the (g-)Askey scheme of 
hypergeometric orthogonal polynomials having discrete orthogonality measures. They are 
also called orthogonal polynomials of a discrete variable [U [6]. For example, they are the 



(g-)Racah, the (g-) (dual) Hahn, the (g-)Krawtchouk, the (g-)Charlier and the (g-)Meixner 
polynomials [H EJ [7] . Various expressions of the transition probability are given explicitly 
together with the totality of the eigenvalues and the measures of the Karlin-McGregor type 
representation [8]. 

It is well-known that the one- dimensional Fokker-Planck equation is related by a similar- 
ity transformation to a corresponding one-dimensional time-independent Schrodinger equa- 
tion [1] or the eigenvalue problem for a suitable Hamiltonian. In other words, solutions of an 
exactly solvable Schrodinger equation give the solutions of the corresponding Fokker-Planck 
equation, which is now exactly solvable. The exact solvability means that the totality of 
the eigenvalues (in these cases, all are discrete) and the corresponding eigenfunctions are 
obtained exactly. Here the Hamiltonian in quantum mechanics is an hermitian (self-adjoint) 
linear operator in a certain Hilbert space. A natural discretisation of the Hamiltonians of 1-d 
quantum mechanics is hermitian matrices of a finite or infinite dimensions. Recently exactly 
solvable 'matrix' quantum mechanics was proposed by Odake and the present author [9] 
by adopting special types of tri-diagonal Jacobi matrices of finite or infinite dimensions as 
Hamiltonians. The eigenfunctions are spanned by the above mentioned orthogonal polyno- 
mials of a discrete variable. The corresponding discretisation of the Fokker-Planck equation 
is, as expected, the birth and death process with a reflecting wa//(s) fl3.20p . Among the 
18 exactly solvable birth and death processes to be explored in this paper, some are quite 
well-known having the linear [21 [31 [10l [11] and quadratic [12] birth and death rates, corre- 
sponding to the Meixner §4.11[ Charlier §4.12[ Krawtchouk §4.41 and Hahn §4.21 polynomials. 
The others have rational functions (of the population x) of the birth and death rates cor- 
responding to the dual Hahn §4.31 and Racah §4.11 polynomials and some others have q ±x - 
linear, quadratic and rational birth and death rates. The most generic one is the g-Racah 
polynomial §4.51 having q x rational birth and death rates f)4.2ip . 

This paper is organised as follows. In section two, the general properties of the Hamil- 
tonians in 1-d quantum mechanics (and/or the hermitian matrices) are reviewed in §2.11 
The relationship between the Schrodinger equation and the corresponding Fokker-Planck 
equation is recapitulated in §2.21 and the solutions of the initial value problem of the Fokker- 
Planck equations and the transition probabilities are expressed in terms of the orthogonal 
polynomials constituting the eigenfunctions of the corresponding Scrodinger equation. In 
section three the birth and death operator is derived from the generic form of the Hamilto- 
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nian of the exactly solvable 'matrix' quantum mechanics of [9]. The solutions of the initial 
value problem of the birth and death equations and the transition probabilities are expressed 
in terms of the orthogonal polynomials constituting the eigenf unctions of the correspond- 
ing Scrodinger equation of the 'matrix' quantum mechanics. Various equivalent expressions 
of the transition probabilities are derived in terms of the dual polynomials. Section four 
provides various data, the birth and death rates, the energy spectra and the sinusoidal coor- 
dinates, the stationary probability, the normalisation constants, and the eigenpolynomials, 
of the exactly solvable 18 models, which are sufficient to evaluate the transition probability 
explicitly. These 18 models are named after the eigenpolynomials, such as the (g-)Racah, 
etc. The final section is for a brief summary and comments. Appendix A provides the col- 
lection of the definitions of basic symbols and functions for self-containedness. Throughout 
this paper we use the parameter q in the range < q < 1. 

2 Fokker-Planck Operator from Hamiltonian 

Here we recapitulate the well-known connection between the Fokker-Planck equation and 
the Schrodinger equation p] in order to introduce appropriate notation and settings for the 
main purpose of the paper; connecting the birth and death process to the 'matrix' quantum 
mechanics to be explored in the next section. 

2.1 Properties of Hamiltonians 

Throughout this paper we discuss one degree of freedom systems only. The Hamiltonians to 
be discussed in this paper are time independent and share the properties listed below. Most 
properties are common to the the Hamiltonians having the continuous dynamical variable 
x (to be used for the Fokker-Planck equation) and the discrete dynamical variable x (to 
be applied to the birth and death processes). They are expressed by the same symbols. 
When they need different symbols, like the L 2 and £ 2 norms, two different expressions are 
shown in a curly bracket as in (12. 3p and (12. 7\i . The upper (lower) one is for the continuous 
(discrete) dynamical variable case. The former (the continuous variable) case corresponds 
to the ordinary quantum mechanics and the 'discrete' quantum mechanics with the pure 
imaginary shifts [13J, which gives rise to the 'deformed' Fokker-Planck equations [14] . 
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(i) Factoris ability, 

H = A^A, (2.1) 

in which t denotes the hermitian conjugation with respect to the standard L 2 (£ 2 ) inner 
product, see ( 12. 3D . This also means that the Hamiltonian 7i is positive semi- definite. 

(ii) Completeness of its eigenf unctions <f> n (x) belonging to discrete eigenvalues (all distinct), 

H(f> n {x) = £{n)(f> n {x), £{0) <8{1) < ■■■ , (2.2) 
and all the eigenvectors are square normalisable and orthogonal with each other 

def | f <f> n (x)*<j> m (x)dx \ 1 
22 x 4>n( X ) 4>m{X) J d n 

The range of the integration (summation) depends on the specific Hamiltonian. Any 
element in the Hilbert space H is expanded by {4>n}'- 

V/eH4/-^W„, n = f cW> n , f n = (^ n ,f). (2.4) 

n 

Here and hereafter / denotes a normalised vector / = // a/ (/, /). We choose all the 
eigenfunctions {0„} to be real, which is always possible in one- dimensional quantum 
mechanics. 

(iii) The groundstate wavefunction <po is annihilated by A and is positive everywhere, 

A(f> (x) = HM%) = 0, S(0) = 0, o (x) > 0. (2.5) 

(iv) The eigenfunction <p n (x) is <f)o(x) times a polynomial, 

<j> n {x) = o (x)P„(r/(x)), n = 0,l,2,..., P = 1, (2.6) 

in which a real function ^(x) is called a sinusoidal coordinate [T5l l9l [T3] . In other words 
P n {v) is an orthogonal polynomial with the orthogonality measure 0o( a; ) 2 

f(j) Q (x) 2 P n (r](x))P m (r](x))dx 1 
E,0o(x) 2 P„(r / (a;))P m (r ? (a;)) J ~ Bm " 

(v) TTie similarity transformed Hamiltonian 7i with respect to 0oOz) 

ft d = (P^oHo O (2.8) 

provides a differential (difference) equation governing the polynomial P n (n(x)). 



d_ 
di 



2.2 Fokker- Planck Equation 

The Fokker-Planck equation in one dimension reads 

: V(x;t) = L FP V(x;t), V(x;t)>0, J V(x; t)dx = 1, (2.9) 

in which V(x; t) is the probability distribution over certain continuous range of the parameter 
x; for example (—00,00), (0, 00) or (0, 7r). The Fokker-Planck operator Lpp corresponding 
to the Hamiltonian H (E21) is defined by [H, [H] 

L FP = -fooHofa 1 , (2.10) 

in which O is defined in ( 12. 51) 1*1 This guarantees that the eigenvalues of L FP are negative 
semi-definite. The square normalised groundstate eigenfunction 4>q(x) provides the stationary 
distribution 4>o{x) 2 of the corresponding Fokker-Planck operator: 

-0 o (x) 2 = L FP Mx) 2 = 0, / ^{xfdx = 1. (2.11) 

It is obvious that (fto(x)<j) n (x) is the eigenvector of the Fokker-Planck operator Lpp'- 

L FP (f> (x)(f) n (x) = -S (ri)(j)o{x)(f) n (x), n = 0, 1,.... (2.12) 

Corresponding to an arbitrary initial probability distribution V{x\ 0), (with J V(x;0)dx = 1), 
which can be expressed as a linear combination of {<j) (x)(t) n (x)}, n = 0, 1, . . ., 



V(x;0) = 4> (x)^2c n (f) n (x), c = 1, c n = (0 n , 4> (x) 1- P(x;0)), n = l,2,..., (2.13) 



def 

C n (p n {X), C = 1, C n 
n=0 

we obtain the solution of the Fokker-Planck equation 



V(x;t) = Mx)Y,Cne- £{n)t 4> n {x), t > 0. (2.14) 

n=0 

This is a consequence of the completeness of the eigen] 'unctions {(p n (x)} (the polynomials) 
of the Hamiltonian 7i. The positivity of the spectrum S (n) > 0, n > 1 (12.21) guarantees that 
the stationary distribution 4>l(x) is achieved at future infinity: 

lim V(x;t) = 4%(x). (2.15) 

t^oo 



1 It should be emphasised that the inverse similarity transformation in terms of </>o is used here: L F p — 



The transition probability from y at t = (i.e., V{x\ 0) = <5(x — y)) to x at t is given by 

oo 

P{y, x; t) = M^Mvr 1 Yl ^ £(n)t Mx)k(y), t > 0. (2.16) 

n=0 

In terms of the polynomial P n (r](x)), (12. 6p . it is expressed as 

oo 

V(y, x; t) = o (x) 2 £ e-Wp^)^^)), f > , (2.17) 

n=0 

in which d n is the normalisation constant (12.31) . (12.71) . 

As shown in p3] in some detail, various examples of exactly solvable quantum mechanics 
[TBI fT7] and the 'discrete' quantum mechanics with the pure imaginary shifts [TH [18j [15] 
provide many explicit cases in which the transition probability (I2.16I) - (I2.17I) can be obtained 
exactly. The corresponding orthogonal polynomials are the Hermite, Laguerre and Jacobi 
polynomials in the ordinary quantum mechanics [TBI [T7] and the Meixner-Pollaczek, con- 
tinuous (dual) Hahn, Wilson and Askey- Wilson polynomials [El [13] and their degenerate 
polynomials, like the continuous g-Hermite polynomials. 

3 Birth and Death process from 'Matrix' Quantum 
Mechanics 

The birth and death equation is a discretisation of the Fokker-Planck equation in one di- 
mension (12.91) . It reads 

^-V(x;t) = (L BD V)(x;t), V(x;t)>0, £)P(a;;t) = 1, (3.1) 

X 

in which V(x; t) is the probability distribution over a certain discrete set of the parameter 
x. Here we simply take a set of consecutive non-negative integers, either finite or infinite: 

x e Z, x G [0, N] or [0, oo). (3.2) 

The exactly solvable birth and death operator or a matrix Lbd is derived from the generic 
form of an exactly solvable Hamiltonian TC of a 'discrete' quantum mechanics with real shifts 

H = -y/B(x) e 9 ^D(x) - ^D(x) e~ 9 y^B(x) + B(x) + D(x), (3.3) 

in which the two functions B(x) and D(x) are real and positive but vanish at the boundary: 

B(x) > 0, D(x) > 0, D(0) = ; B(N) = for the finite case. (3.4) 
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The explicit forms of the functions B(x) and D(x) are given in each subsection of section 
four, which are named after the orthogonal polynomials appearing as the main part of the 
eigenfunctions. In the Hamiltonian (13.31) e ±d are formal shift operators acting on a function 
/ of x as 

(e ±9 f)(x) = f(x±l). 

Thus the Schrodinger equation 7iijj(x) = £ip(x) is a difference equation with real shifts: 



(B(x) + D{x))i){x) - ^B{x)D{x + l)i){x + 1) - y/B(x - l)D(x) ip{x - 1) = Sip{x), 

x = 0,l,...,(N),.... (3.5) 

The boundary condition D(0) = is necessary for the term ip{— 1) not to appear, and 
B(N) = is necessary for the term ip(N + 1) not to appear in the finite dimensional matrix 
case. 

Although the Hamiltonian 7i (13.31) is presented in a difference operator form, it is in fact 
a real symmetric tri-diagonal (Jacobi) matrix: 

T~L = C%e,i/)) 7~tx,y = 7~Ly,x, (3.6) 

H x , y = -y/B(x)D(x + l)6 x+1 , y - y/B(x-l)D(x)5 x - lt y + (B(x) + D(x))6 x>y . (3.7) 

As mentioned above, the Hamiltonian is factorisable (12. ip . TC = A^A: 

A^ = y/B(x) - y/D{x) e~ 9 , A= ^fB{x) - e d ^fD{x). (3.8) 

In the matrix form A^ has the diagonal and sub-diagonal elements only and A has the 
diagonal and super-diagonal elements only 

{A v ) x , y = a/ B{x) 8 x<y - a/ D(x) 5 x _ lt y, A x>y = \/ B(x) 5 x>y - y/ D(x + 1) 6 x+1>y . (3.9) 

The equation (12.51) determining the groundstate wavefunction <J)q is easy to solve, since 
.4.00 = is a two term recurrence relation: 



Ux + l) B(x) 



0o(x) y D(x + 1) 
It can be solved elementarily with the boundary (initial) condition 0o(O) = 1, 



(3.10) 



00 0*0 
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With the standard convention Yik=n * = 1> the expression (13.111) is valid for x = 0, too. For 
the infinite matrix case, the requirement of the finite £ 2 norm of the eigenvectors 

OO OO X— 1 „/ \ 

z=0 z=0 y=0 Vy ' 

imposes constraints on the asymptotic behaviours of B(x) and D(x). 

With the above explicit form of the groundstate wavefunction ^o^); the similarity trans- 
formed Hamiltonian (I2.8P is easily obtained 

H = <f) 1 o H o O = B(x)(l - e a ) + D(x)(l - e~ d ). (3.13) 

As mentioned above, Ti provides the difference equation for the polynomial eigenfunctions, 

(HP n )( V (x)) = S{n)P n ( V (x)), (3.14) 

that is, 

B(x) (P n (v(x)) - P n {r){x + 1))) + D(x) (P n ( V (x)) - P n ( V (x - 1))) = £(n)P n ( V (x)). (3.15) 

The eigenpolynomials {P n } are the orthogonal polynomials of a discrete variable. See §5 
of [H] for various forms of B(x) and D(x) and the corresponding orthogonal polyno mails. 
It is also recapitulated in section 4 of this paper. For example, they are the (g-)Racah, 
the (g-)( dual) Halm, the (g-)Krawtchouk, the (g-)Charlier and the (g-)Meixner polynomials 
[H El U\ • As a matrix, Ti is another tri-diagonal matrix 

H = (H Xt y), H x ,y = B(x)(5 X} y - S x+ i tV ) + D(x)(S x>y - S x -. 1>y ) . (3.16) 

Corresponding to (I2.10p . the inverse similarity transformation of the Hamiltonian TC 
supplies the birth and death operator Lbd'- 

Lbd = -0o o H o O 1 = (e~ d - l)B(x) + (e 9 - l)D(x). (3.17) 

Obviously the stationary probability is given by (po(x) 2 = dl4> (x) 2 . In the matrix form, L B d 
is again tri-diagonal: 

L B d = (L B D x ,y), L B D x , y = B(x - - B(x)5 X) y + D(x + l)5 x+liy - D{x)8 x>y . (3.18) 

In fact, —Lbd is the transposed matrix of Ti: 

~ L B d = (T~t) , —L B d x , v = 7~L y>x . (3.19) 
8 



With the explicit form of the birth and death operator Lbd, the birth and death equation 
(13. ip in our notation reads 

d 

—V{x;t) = ^2L BDXiy V(y;t) 
v 

= ~{B(x) + D(x))V(x; t) + B(x - l)V{x - 1; t) + D(x + l)V(x + 1; t). (3.20) 

The standard interpretation is that x is the population of a group, V(x; t) is the probability 
for the group to have the population x at the time t, and B(x) is the birth rate, D(x) is the 
death rate, respectively, when the population is x. It is quite easy to remember. This is to 
be compared with the standard notation, for example, [8], XVII. 5 of [2], §5.2 of [5]: 

d 

Q^Pnyt) = _ (An + H>n)Pn[t) + X n -lPn-l(t) + fJ, n +lPn+l (t) , (3.21) 

in which A n is the birth rate and /i„ is the death rate. The following translation table of the 
notation will be helpful. 





standard [21 [5] 


this paper 


population 


n = 0,l,...,(iV),... 


x = 0,l,...,(iV),... 


probability 


Pn(t) 


V(x;t) 


Birth rate 


A„ (Ajv = 0) 


B(x) {B(N) = 0) 


Death rate 


(yUo = 0) 


D{x) (D(0) = 0) 



Table I: Translation Table. 



The boundary condition for the finite case, A at = (B(N) = 0) ( 13.41) is said that the system 
has a reflecting wall at the population N. 

The transition probability from y at t = (i.e., V(x; 0) = 8 x>y ) to x at t has exactly the 
same expression as that in the Fokker-Planck equation (12.161) 

V(y, x- 1) = Mx^y)- 1 ^ e- s ^ n {x)^ n {y), t > 0. (3.22) 

n=0 

In terms of the polynomial P n (r)(x)), (\2.6L it is expressed as 

V(y, x; t) = M*) 2 E d n e- £(n)t P n (v(x))Pn(v(y)), t > 0. (3.23) 

n=0 

It should be emphasised that in these formulas (13 .22[) - (I3.23j) everything is known including 
the measure in contradistinction to the general formula by Karlin-McGregor [8]. 
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Let us mention several equivalent expressions of the transition probability (I3.23j) in terms 
of the dual polynomials [19], [201 12H [9] . It is well-known that with proper normalisation 

^(0) = = £(0), P = 1 = Q , P n (0) = g,(0) = 1, (3.24) 

the two polynomials, {P n (r])} and its dual polynomial {Q x (£)}, coincide at the integer lattice 
points [9]: 

P n (r](x))=Q x (£(n)), n = 0, 1, . . . , (N), . . . , x = 0, 1, . . . , (N), . . . . (3.25) 

The dual polynomial {Q x (£(n))}, x = 0,1, . . ., is a right eigenvector of the similarity trans- 
formed Hamiltonian 7i matrix with the eigenvalue £{n): 

n x , y Q y (£(n)) = £{n)Q x (£(n)). (3.26) 

y 

The above equation is the three term recurrence relation for the dual polynomials {Q x (£)}: 

(B(x) + D(x)) Q x (£(n)) - B(x)Q x+1 (£(n)) - D(x)Q x ^(£(n)) = £(n)Q x (£(n)), (3.27) 

Q = 1, Q x {£) = (B(0) - £)/B{0), Q 2 {£) = (B(0) - £)(B(1) + D(l) - £)/{B{0)B{l)), .... 

(3.28) 

For historical reasons, this polynomial Q x (£) is called the birth and death polynomial or the 
Karlin-McGregor polynomial [8]. 

In terms of the dual polynomials or the Karlin-McGregor polynomial, the transition 
probability is 

V(y, x; t) = o (x) 2 ^ d 2 n e- e WQ x (E(n))Q y {S{n)), t > 0. (3.29) 

71=0 

Following [5], let us introduce 

F x {£{n))^ <p {x) 2 Q x {£{n)). (3.30) 

Since L BD and TC is related by 

Led = -0o ° H ° <Po 2 , (3-31) 

it is easy to see that F x (£(n)) is a left eigenvector of TC and thus a right eigenvector of the 
birth and death operator Lbd'- 

y^ j L B Dx,yFv(E{n)) = -(j)o{xf^H x , y Q y {£{n)) 
y y 

10 



= -£{n)<p {x) 2 Q x {£{n)) = -£{n)F x {£{n)). (3.32) 



In terms of the right eigenvectors of Lbd-, we obtain another expression of the transition 
probability |5j 

V(y,x-t) = -^—Y,d 2 n e- £ ^ t F x (£(n))F y (£(n)), t > 0. (3.33) 

My) 2 fo 

The explicit forms of the transition probability (I3.22p . (13.231) . f!3.29j) and f!3.33|) can be 
evaluated straightforwardly if the Hamiltonian 7i of an exactly solvable discrete quantum 
mechanics is given. Thus we may call the functions B(x) and D(x) in the Hamiltonian 7i of 
an exactly solvable discrete quantum mechanics (13.31) . the birth and death rates of an exactly 
solvable birth and death process. As mentioned above, the association of the birth and death 
rates and the orthogonal polynomial in this paper and in the literature jH [12] are dual to 
each other. Therefore the names of the polynomials in the next section are the dual of the 
corresponding Karlin-McGregor polynomial except for the self-dual cases of the Krawtchouk 
§13 Meixner gH and Charlier 

In the subsequent section we will present 18 examples of exactly solvable birth and death 
processes. 

4 18 Examples 

Now let us proceed to give the 18 explicit examples of exactly solvable birth and death 
processes. The input is simply the function forms of the birth and death rates B(x) and 
D(x). The rest is calculable. But here we also provide other data, taken from [9], such 
as the energy eigenvalue £(n), the sinusoidal coordinate r)(x), the unnormalised stationary 
probability <fio(x) 2 , the normalisation constants d 2 n and the polynomials P n {v)- Following 
the order of our previous work on the exactly solvable discrete quantum mechanics [9J, we 
handle the most generic one first, and then followed by the simpler ones. There is a logical 
reason for this order. The simpler ones are usually obtained by specialising or restricting the 
parameters of the generic ones. Each example is called by the name of the corresponding 
orthogonal polynomial P n {v) with the number e.g. [KS3.2] attached to it indicating the 
subsection in the standard review of Koekoek and Swarttouw [7]. The finite (N) cases are 
discussed first and then the infinite ones. In each group the Askey-scheme of hypergeometric 
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orthogonal polynomials (non-g polynomials) will be discussed first and followed by the q- 
scheme polynomials. 

Please note that the set of parameters is slightly different from the conventional ones 
[HO [7] for some polynomials, the reason explained in For some polynomials, for example, 
the (q-) Racah, (dual, q-) Hahn, etc, there are many non-equivalent parametrisations of B (x) 
and D(x), which could lead to non-equivalent birth and death processes. Here we give only 
one of them as a representative, since the purpose of the paper is to show exactly solvable 
structure, not to provide an exhaustive list of all solvable models. See [9] for more general 
parametrisations and the allowed ranges of the parameters. In the same spirit we did not 
include some of the polynomials listed in [9]. 

Finite Dimensional Cases 

4.1 Racah [KS1.2] 

The Racah polynomial is the most generic hypergeometric orthogonal polynomial of a dis- 
crete variable. All the other (non-g) polynomials are obtained by restriction or limiting 
procedure. The function B(x) and D(x) depend on four real parameters a, b, c and d, with 
one of them, say c, being related to N, c = —N: 

. . (x + a) (x + b) (x + c) (x + d) (x + d — a)(x + d — b)(x + d — c)x 

{X} = (2x + d)(2x + 1 + d) ' {X) = (2x-l + d)(2x + d) ' 

(4.1) 

The other data are: 

S(n) = n(n + d), r](x) = x(x + d), d= a + b + c — d — 1, (4-2) 

a>b, d>0, a> N + d, < b < 1 + d, (4.3) 

Ux)^ (Q ' & ' cA (4 4) 

rov ; (l + d-a,l + d-6,l + d-c,l) a . d ' y ' ' 

2 (a,b,c,d) n 2n + d (—l) N (l + d — a,l + d — b,l + d — c)N 
d n = = — x . 

(1 + d— a, 1 + d — b, 1 + d — c, l) n d (d + l) N (d + 1) 2 n 

(4.5) 

Here (a) n is the Pochhammer symbol flA.lj) . Throughout this section, the format for d\ 
consists of two parts separated by a x symbol: d\ = (d^/d^) x d^. The second part d^ 
satisfies the relation J2 x 4>o( x ) 2 — l/^o- The polynomial is 

Pn(v(x))=M~ n,n + \~ X,X + d X )' ( 46 ) 

V a, b, c J 
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in which 4F3 is the standard hypergeometric series (1A.3j) . The dual polynomial is again the 
Racah polynomial with the parameter correspondence (a, b, c, d) <-> (a, b, c, d). The rational 
(a quartic polynomial divided by a quadratic polynomial) birth and death rates (14. II) have 
not yet been discussed but the Racah polynomial appears in [12]. 

4.2 Hahn [KS1.5] 

This is a well-known example of quadratic (in x) birth and death rates with two real positive 
parameters a and b: 

B(x) = (x + a)(N-x), D(x) = x{b + N-x). (4.7) 

It has a quadratic energy spectrum 

E(n)=n(n + a + b-l), r,(x) = x, o (x) 2 = g , ^t|^ , (4.8) 

2 = N\ {a) n {2n + a + b-l){a + b) N (b) N 

n n\(N-n)\ (b) n (n + a + b-l) N+l X (a + b) N ' l " J 

. , . . /—n,n + a + b—l,—x \ 

P n ( V (x)) = 3 F 2 ( Ar 1. (4.10) 



a, -N 

The dual polynomial is the dual Hahn polynomial of the next subsection 14.31 The quadratic 
birth and death rates are discussed in [T2] associated with the dual Hahn polynomial. 

4.3 dual Hahn [KS1.6] 

The set of parameters is the same as the Hahn polynomial case. The birth and death rates 
are rational functions of x, 

, {x + a){x + a + b- 1){N - x) _. , x(x + b - l)(x + a + b + N - 1) 

B(x) = —, TT7 tt^. D(x) = - 1 , w TV, (4.11) 

v ; (2x-l + a + 6)(2x + a + 6) ' v 7 (2x - 2 + a + b)(2x - 1 + a + b) ' v 7 
giving rise to a linear energy spectrum 

BV-n. rj(x)=x(x + a + b-l), -^fcJ^^Z^, 

(4.12) 

,2 = M (a)n 0)jV-n (&)jV 1 q\ 

% n\{N-n)\ (b) N (a + b) N > { ' 

„ . , ., „ /— n, x + a + 6 — 1, — x \ 

P n ( V (x)) = 3 F 2 ( a _iv V- ( 4 - 14 ) 
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4.4 Krawtchouk [KS1.10] (self-dual) 

The case of linear birth and death rates are a very well-known example (the Ehrenfest model) 
[TT] of an exactly solvable birth and death processes [2j [3]: 

B(x)=p(N-x), D(x) = (l-p)x, 0<p<l, (4.15) 

£{n) = n, rj(x) = x, (4-16) 

= (JL.)-, 4= *"(^_)- „(!_„)» (4.17) 

x! (A* — x)! VI — pJ n\{JS — nj! VI — pJ 

Pn(v(x)) = 2F 1 (~^~ X \p- 1 ). (4.18) 

This is a simplest example of self-dual polynomials. The stationary probability <f)o(x) 2 dl is 
the binomial distribution. 

4.5 g-Racah [KS3.2] 

This is the first example of the g-scheme of the orthogonal polynomials. Among them the q- 
Racah polynomial is the most generic. The set of parameters is four real numbers (a, b, c, d), 
which is different from the standard one in the same manner as for the Racah polynomial. 
We restrict them 

c = q~ N , a<b, 0<d<l, 0<a<q N d, qd < b < 1, d<q~ l , d = abed" 1 q~ l . (4.19) 

The functions B(x) and D(x) are 

p/ \ _ (1 ~ «<f X 1 - bq x )(l - cq x )(l - dq x ) 
[X) ~ (l-dq 2 *)(l-dq 2 *+l) ' 1 J 

= j (1 ~ a-^tf )(1 - fc^Kl - c^Kl - q x ) 
1 J ' (1 -dq 2 *- 1 )^ -dq 2x ) ' 1 ' ] 

The other data are 

S(n) = (q~ n - 1)(1 - rfg"), r/(x) = (g~* - 1)(1 - dg»), (4.22) 
,,,2 (a,b } c,d;q) x 1 - dq 2x 

(po{x) = = — , (4.23) 

(a~ 1 dq,b~ 1 dq,c~ 1 dq,q;q) x d x i — d 

j2 (a,b,c,d;q) n 1 - dq 2n {-l) N {a~ 1 dq,b- 1 dq,c~ 1 dq;q) N d N q l 2 N ( N +V 

X 



(a 1 dq, b 1 dq,c 1 dq,q;q) n d n l-d (dq;q) N (dq;q) 2 N 

(4.24) 
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Pn(v(x))= 4 <j> 3 ( q ~ n ' dqn ' q ~ X ' dqX g;g), (4.25) 



a, b, c 

in which 4 3 is the basic hypergeometric series ()A.4j) and (a; g) n is the g-Pochhammer symbol 



( 1A.2j) . The dual g-Racah polynomial is again the g-Racah polynomial with the parameter 
correspondence (a, b, c, d) <-> (a, 6, c, d). 

4.6 g-Hahn [KS3.6] 

The g-Hahn polynomial has two positive parameters a and b and the birth and death rates 
are quadratic polynomials in q x : 

B{x) = {I - aq x ){q x ~ N - 1), D(ar) = ag- x (l - q x ){q x ' N - b), < a,b < 1. (4.26) 

The other data are 

£(n) = {q- n -l){l-abq n - 1 ), 77(a) = q~ x - 1, (4.27) 

0o (x) = t r7 r 77 — v , 4.28 

{q;q)x{q;q)N-x {b;q) N a x 

d 2 = (Q',q)n (a,abq- 1 ;q) n 1 - afeg 2 "' 1 ^ (b;q) N a N 

{q;q) n {q;q)N-n {abq N ,b;q) n a n 1 - abq- 1 {ab;q) N 

PM*)) = 302 ( 9 ~"' a6< C' 9 ~* ?;<?)■ (4-30) 

Obviously the g-Hahn and dual g-Hahn are dual to each other. 

4.7 dual g-Hahn [KS3.7] 

For obvious reasons, we adopt the same parameters (a, b) for the g-Hahn and dual g-Hahn 
polynomials. The birth and death rates are rational functions of q x : 



£(n)=g- n -l, r/(x) = (g^ - 1)(1 - afeg^ 1 ), (4.33) 

<p (x) 2 = (iigk (a, abq' 1 ;q) x 1 - a&g^" 1 
'" {q;q) x {q;q) N - x (abq N ,b;q) x a x I - abq~ l 
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d 



(q;q) N (a;q) n (b;q) N ~n (b;q) N a 

X 



N 



(q;q)n(q;q)N-n {b;q) N a n 

( q n ohq x ~^ q x \ 
Pni'nix)) = 30 2 f ' _ JV ' q;qj. 

4.8 quantum g-Krawtchouk [KS3.14] 



(ab; q) 



N 



(4.35) 
(4.36) 



This has one positive parameter p > q N . The birth and death rates are quadratic polyno- 
mials in q x : 



B(x) = p~W~ N - 1), D(x) = (1 - <f )(1 - p^q^' 1 ), 
£(n) = 1 - q n , r)(x) = q~ x - 1, 

2 (Q',Q)n p-XqX(x-l-N) 



dl 



(q ; q)x(q ; q)N- x (p- 1 q~ N ; q) x ' 
(g;g)jy 

(q;q)n(q;q)N-n (p~ 1 q~ n ;q) n 



p n q Nn , _i _at N 



q n , q x 

n -N 



q-pq 



(4.37) 
(4.38) 

(4.39) 
(4.40) 
(4.41) 



4.9 g-Krawtchouk [KS3.15] 

This has one positive parameter p > and the birth and death rates are linear in q x : 

B(x) = q x ~ N - 1, L>(x) =p(l -g x ), 
5(n) = - 1)(1 77(2;) = g"* - 1, 

o (x) 2 



(g j gW r) -21 n \x{x-l)-xN 

I q 1 



d 



PnHx)) 



(q;q)x(q;q)N~x 
2 (q ; 7)v 



(-p;?)t 



1 + pg 2n p N q \ N ( N + l ) 

X 



(q;q)n(q;q)N-n (-pq N + 1 ;q) n p n q^ n( - n+ v l+p ^ 



3^2 



q n , q x , —pq n 

-N 



q-»,0 



(4.42) 
(4.43) 
(4.44) 

(4.45) 
(4.46) 



4.10 affine g-Krawtchouk [KS3.16] (self-dual) 

This has one positive parameter p and the birth and death rates are quadratic polynomials 



m q 



B(x) = {q x ~ N -l)(l-pq x+1 ), D{x) =pq x - N {l-q x ), S < p < q~\ (4.47) 
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£{n) = q n - 1, T](x) = q x ■ 

2 (9 ; q)n ipq ; ?)* 



00 



(q;q) x (q;q) N -x ipq) 3 



N 



(gg ; g)r 

(q;q)n(q;q)N-n ipq) n 



(9; 9) 



Pn{f]{x)) = 302 



q n , q x , 



pq, q 



-N 



q;q 



(4.48) 

x (4.49) 



(4.50) 



Infinite Dimensional Cases 

In contrast to the finite dimensional case, the structure of the polynomials is severely con- 
strained by the asymptotic forms of the functions B(x) and D(x) (13.121) . 

4.11 Meixner [KS1.9] (self-dual) 

This is the best known example of exactly solvable birth and death processes [10] and the 

birth and death rates are both linear in x with simple linear energy spectra Sin) = n and 
r]{x) = x. It has two positive parameters (3 and c: 

B(x) = ^—(x + (3), D(x) = —^—x, (3>0, 0<c<l, (4.51) 
1 — c 1 — c 

£{n) = n, r](x) = x, (4-52) 

0o(x) 2 = ^, ^ = x (1 - cf, (4.53) 
x\ n\ 

P n ( V (x)) = 2 F 1 (~ n, ~ X 1-c- 1 ). (4.54) 



4.12 Charlier [KS1.12] (self-dual) 

This is another best known example of exactly solvable birth and death processes with a 
constant birth rates a > and a linear death rates: 

B(x) = a, D(x) = x, (4.55) 
S{n) = n, r](x) = x, (4.56) 

0o(x) 2 = ^, dl = -^e~\ (4.57) 
x\ n\ 

P n (r ? (x) = 2 F f" n ' ~ X -a' 1 ). (4.58) 



The stationary probability <fio(x) 2 dl (I4.57P is the Poisson distribution. 
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4.13 little g-Jacobi [KS3.12] 

This has two parameters a and b. The birth and death rates grow exponentially as x tends 
to infinity: 

B(x) = a(q- x -bq), D(x) = q~ x - 1, < a < q~\ b<q~\ (4.59) 
E(n) = (q~ n - 1)(1 -abq n+1 ), v {x) = l-q x , (4.60) 

MX)2 = J^t {aq)X > (461) 
d2 = {bq, abq ; q) n a n q n2 1 - abq 2n+l ^ (aq;q)oo ^ ^ 

{q,aq;q) n 1 - abq {abq 2 -^)^' 

p n ( v (x)) = (-a)-v^ ( " +1) ^^ 2 0i( g " n ' ahqn+1 g;g x+1 ). (4.63) 

{bq;q) n \ aq ) 
The normalisation of the polynomial is different from the conventional one. 

4.14 g-Meixner [KS3.13] 

This has two positive parameters b and c. The birth and death rates are quadratic in q x and 
as x goes to infinity, the birth rates tend to zero and the death rates tend to unity: 

B{x) = cq x {\ -bq x+1 ), D(x) = (1 - <f)(l + bcq x ), 0<b<q-\ c> 0, (4.64) 

£(n) = l-q n , r](x)=q- x -l, (4.65) 

Mxf = (frig)* f fr*-* d , = ibq-q) x (-bcq^ 
{q,-bcq;q) x (q, —c~ Y q ; q) n (-c;g)oo 



PnWx)) = 2 0i ( q ' n bq q X | Q ; -c~y +1 ) . (4.67) 
4.15 little g-Laguerre/Wall [KS3.20] 

This has one positive parameter a and both the birth and death rates grow exponentially as 
x tends to infinity: 

B(x) = aq~ x , D(x) = q~ x - 1, < a < q~\ (4.68) 

E(n) = q- n - 1, V {x) = l-q x , (4.69) 

^ = & <=(iS^ X ^-' (4 - 70) 



Pn(i?(x)) = 2 0o( g ~V 3 | 9; a "V)- (4.71) 
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The normalisation of the polynomial is different from the conventional one. 



4.16 Al-Salam-Carlitz II [KS3.25] 

This has one positive parameter a and the birth and death rates are quadratic in q x . As x 
goes to infinity the birth rates tend to zero and death rates tend to unity: 

B(x) = aq 2x+ \ D(x) = (1 - q x )(l - aq x ), < a < q~\ (4.72) 

£(n) = l-q n , r)(x)=q- x -l, (4.73) 



M*f = , ^ , , ^ = p^-x(aq;q) 0O , (4.74) 
P n (v(x)) = 2<t>o( q X 9;a"V)- (4-75) 



The normalisation of the polynomial is different from the conventional one. 

4.17 alternative g-Charlier [KS3.22] 

This has one positive parameter a. The birth rates are constant a whereas the death rates 
grow exponentially as x goes to infinity: 

B(x) = a, D(x) = q~ x - 1, a > 0, (4.76) 

£(n) = (q~ n -l)(l + aq n ), V (x) = l-q x , (4.77) 
^( x y = a x q^ x+1) d2 = a"g^(3"-i) (_ a ^ 1 + aq ^ 1 

(q;q) x ' n {q;q)n (-ag n ;g)oo l + a (-a? ;?)<»' 

Pn(v(x)) = q nX 2^{ q ~ n Q q ' X I q;-a^q- n+1 ). (4.79) 
The normalisation of the polynomial is different from the conventional one. 

4.18 g-Charlier [KS3.23] 

This has one positive parameter a and as x goes to infinity the birth rates tend to zero and 
the death rates tend to unity: 

B{x) = aq x , D(x) = l-<f, a > 0, (4.80) 

E(n) = l-q n , r]{x) = q~ x - 1, (4.81) 

(q;q) x ' dn {-a~ 1 q,q;q) n X (-a;g)«> ' ( ' 4 ' 82 ' ) 

19 



P n (r)(x)) = 20i 



q n , q x 



q ; -a~ x q n+1 



(4.83) 



5 Summary and Comments 

Following the simple line of arguments summarised in the following diagram, we presented 
18 models of exactly solvable birth and death processes and their solutions, the transition 
probabilities. In the diagram 'ES' stands for Exactly Solvable. 



ES Id Quantum 
Mechanical systems 



give solutions 



ES Id Fokker-Planck 
equations 



discretisation 



discretisation 



ES 'Matrix' Quantum 
Mechanical systems 



give solutions 



ES Birth and Death 
processes 



The exactly solvable 'matrix' quantum mechanics, or the 1-d 'discrete' quantum mechanics 
with real shifts was explored in detail in [U] to cover most of the hypergeometric orthogonal 
polynomials of a discrete variable in the (q-) Askey scheme jU El [7]. For the 'explanation' 
of the exact solvability, see a recent work [22] ■ By comparing the present simple results 
with those in the literature [SI El El 12] one would realise the essential role played by the 
energy spectrum £ (n) and the sinusoidal coordinate i](x). They are the eigenvalues of the two 
operators, called the Leonard pair, which characterise the orthogonal polynomials completely 

[13 [201 EI]- 

In this paper we did not discuss the generalisation of the birth and death processes which 
has fio > (-D(O) > 0), the non-vanishing death rate at zero population, although this has 
led to a new type of orthogonal polynomials in the cases when the birth and death rates 
B(x) and D(x) are linear and quadratic in x, [23, 24J. It would be interesting to try further 
generalisation in this direction for which B(x) and D(x) are rational, e.g. the Racah case 
§4. II or g-linear, e.g. the g-Krawtchouk §4.91 or g-quadratic, e.g. the the affine g-Krawtchouk 
§4.101 or even the g-rational, e.g. the g-Racah §4.51 cases. 

It is a big challenge to try and find a closed form expression for 



(5.1; 



n=0 
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appearing as a part of the transition probability (12.171) . (13.231) for various examples in section 
four. To the best of our knowledge, such expressions are known only for the linear energy 
spectrum S(n) oc n. For example, for the Fokker-Planck equation corresponding to the 
harmonic oscillator Hamiltonian, or the Ornshtein-Uhlenbeck process [H [H] , we have: 

d 2 



n 

V(y,x;t) 



dcf 



dx 2 

„2 OO 



x 2 - 1, L FP = ^— + £(n)=2n, jj(x) 

dx 2 dx 



H n(x)H n (y) e _ 2nt 



1 



2 n n\ 



-it 



exp 



(x — ye 2t ) 2 



n=o - yWT-e" 
The last equality was derived based on (6.1.13) of [I]. Another example is 



At 



n d = i -S, + x 2 + 9 -^ 



dx 2 



x* 



[l + 2g), L FP = 
S n = An, r)(x) =x 2 , (3 = g- 1/2, 



^ 2 + 2 ^ ( ^ ^ 
dx 2 dx x ' 



n=0 



2x 2 3 



1-e 



-At' 



exp 



T(n + /3 + l) 
(x 2 +y 2 e- 4t ) 



1-e 



-Ant 



(xye 



-2t\-p 



/ 2xye 



-it 



1 



-if 



(5.2) 
(5.3) 



(5.4) 
(5.5) 
(5.6) 

(5.7) 



in which Jg is the modified Bessel function of order (3. The last equality was derived based on 
(6.2.25) of [1]. We would like to ask experts in special functions and orthogonal polynomials 
to derive such bilinear generating functions for various energy spectra: 



£{n) = n(n + d), q 



1, 1 



q\ (g-»-l)(l-dg»). 



(5i 
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Appendix A: Some definitions related to the hypergeo- 
metric and (/-hypergeometric functions 

For self-containedness we collect several definitions related to the (g-)hypergeometric func- 
tions [7]. 
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o Pochhammer symbol (a) n : 

(a)n = f[(a + k - 1) = a{a + 1) • • • (a + n - 1) = ^7^- (A ' X) 

fe=i ^ 

o g-Pochhammer symbol (a;q) n : 

(o 5 ?)n = f[(l - = (1 - a)(l - ag) ■ ■ • (1 - ag"" 1 ). (A.2) 

fc=i 

o hypergeometric series r F s : 

OO 

CL\ , • • • , d. 



h, ■■■ ,b s 



z 



def 1^1 j ' ' ' ; &r)n Z n , . , 

-^(6i,-,6 8 )„^!' (A - 3) 



where (ai, • • • , a r ) n = f ]lj=i( a i)n = ( a i)« • • • (a r )n- 

o g-hypergeometric series (the basic hypergeometric series) r (f> s : 

/dl, • • • , a r \ def >^ (ai, ■■■ ,a r ;g)» \(l+ S -r)n (l+ S -r)n(n-l)/2_£]|_ /a ^ 

r *'W,'--,b. q ' Z )-^(h,...,b.;q) n ( V q (g;g)„' ^ 

where (ai, • • • , a r ; g) n = f rij=i( a j 5 0)n = («i ; q)n • • • iflr 5 ?)n- 
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